Systems with two symmetric absorbing states: relating the microscopic dynamics with 

the macroscopic behavior 
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We propose a general approach to study spin models with two symmetric absorbing states. Start- 
ing from the microscopic dynamics on a square lattice, we derive a Langevin equation for the time 
evolution of the magnetization field, that successfully explains coarsening properties of a wide range 
of nonlinear voter models and systems with intermediate states. We find that the macroscopic 
behavior only depends on the first derivatives of the spin-flip probabilities. Moreover, an analysis 
of the mean-field term reveals the three types of transitions commonly observed in these systems 
-generalized voter, Ising and directed percolation-. Monte Carlo simulations of the spin dynamics 
qualitatively agree with theoretical predictions. 
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Introduction. The non-equilibrium dynamics of inter- 
acting particle systems with absorbing states is a central 
issue in modern Statistical Mechanics [EH- In the last 
years much interest has been given to the special case 
of models with two symmetric (^-symmetry) absorbing 
states [3, 13, [!{, due to the relevance in many different 
contexts of a subclass of them, the so-called voter mod- 
els. They have been widely used in diverse disciplines to 
study different dynamics, such as species competition |6|, 
allele frequency in genetics 0] , kinetics of heterogeneous 
catalysis 8| , and more recently, opinion formation Q and 
language spreading [Icj. In these models, the state of a 
particle in a lattice site evolves according to the density 
of states in its near neighborhood. Interesting dynam- 
ical behaviors arise depending on the specific updating 
rules, the number of states and the functional form of 
the transition probabilities between configurations. For 
instance, it has been found that the addition of memory 
or inertia in the spin dynamics lH E2] , the introduction 
of intermediate states 3, 14, 15 1, or the use of non-linear 
transitions [li} , result in a drastic change of the coarsen- 
ing properties and final outcome of the system. 

Despite that the dynamical rules of the models are very 
different in nature, many of them seem to share the same 
macroscopic behavior, such as coarsening and criticality. 
However, the minimal conditions that a microscopic dy- 
namics must hold in order to observe a particular behav- 
ior have not been clearly identified yet. In other words, 
given a spin model defined by its flipping transition prob- 
ability and interaction range, can we anticipate how the 
system will evolve over time ? 

In this article, we try to answer this question by de- 
veloping an approach that connects the microscopic dy- 
namics with the macroscopic space-time evolution of the 
system in square lattices. We derive a Langevin equation 
for the magnetization field, and find that, at the macro- 
scopic level, the properties are only determined by the 
first three derivatives of the transition probabilities. 

Our Langevin equation coincides with that postulated 



by Hammal et. al. 3] by symmetry arguments, but now 
the coefficients of the different terms have a clear expla- 
nation in terms of the transition probabilities. The anal- 
ysis of this equation helps to understand some of the open 
questions about phase ordering in these systems, that is, 
whether the coarsening is driven by curvature like in the 
Ising model [l7j], or it is without surface tension like in 
the original voter model (VM) jj4j. This approach also 
explains, from a different perspective than in [151 ]. why 
adding intermediate states to the VM leads to an effec- 
tive surface tension. Moreover, numerical simulations of 
the spin dynamics with different interaction ranges con- 
firm the three possible classes of phase transitions un- 
veiled in Q , and clarifies an apparent controversy found 
between previous works 

an. 

The model and the Langevin equation. Each site r = 
(n, .., r c i) of a (i-dimensional square lattice is occupied by 
one particle with a spin that can assume either value 1 
(up) or —1 (down). The dynamics consists of choosing, 
at each time step, a site r at random and flipping the 
spin S r at this site with a probability that is a function, 
f(—S T ip r ), of the product between S r and the particle's 
local magnetization ip r = \ ~}2 r i / r S r i , where the sum is 
over the z nearest-neighboring sites r' of r. z is an arbi- 
trary integer number that defines the interaction range 
(e.g. z = 4 for first-nearest neighbor interactions). In or- 
der to ensure that the fully ordered configurations S r = 1 
or S r = — 1 for all r are absorbing, the flipping probabil- 
ity f must vanish when the spin is aligned with all its 
neighboring spins, i.e, /(— 1) = 0. 

We want to derive a Langevin equation for the field 
4> r (t), that is a continuous representation of the spin 
at site r, at time t. For this we follow a standard ap- 
proach (see HH), and consider an ensemble of Q copies 
of the system, each copy representing a particular spin 
configuration. This is equivalent to assume f2 spin parti- 
cles at each site of the lattice (our microscopic model 
corresponds exactly to Q = 1, but this substitution 
can be made at the end of the calculation). In this 
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approach (f> r (t) is replaced by the average spin value 
4> r (t) — > nS^Li^r) ana ^ by the average local field 
?/> r — ► ^Er'/r^'Wi at s ^ e r at time t. In a time 
step, a site x and one particle from that site with spin 
S x are randomly chosen. Then the spin attempts to 
flip with probability f(—S x ip x ). If the flipping occurs, 
then the field on the entire lattice, represented as {(/)}, 
changes only at site x by — 2S*/£l. Thus the rising and 
lowering transition rates are W ({</>} — > {</>} ± f|^ x . r ) = 
W ± (c)), x, t) = | (1 =F ^x) f(±ijj x ). We can write a master 
equation for the time evolution of the probability distri- 
bution V({4>}, t), which after an expansion to second or- 
der in l/fi leads to the following Fokker-Planck equation 
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(1) 

■<j>,v,t)-w-{4>,v,t)\v{{4>},t)] 

,r,t) +W-(<j>,r,t)]V({<p},t)}. 



From the above equation, using the expressions for 
W^fy, r,t) and rescaling time with f2, we arrive to the 
Langevin equation 



at 



= [1 - Mt)} f(A) - [1 + Mt)} f(-A) + vAt), 

(2) 

where 7] r (t) is a Gaussian white noise with correlations 

+ [1 + r (i)] /(-^ r )}^, r ^(t - i')/r! 1/2 . (3) 

Note that up to now our derivation is completely general, 
and the fact that the system has two absorbing states is 
only present in the condition for the flipping probability. 
Now we look for an approximation to Eq.([2]) which, how- 
ever, captures the behavior of a wide range of absorbing 
Z2 models. As explained in [3j, this is obtained with a 
</> 6 model, i.e. when the right hand side of Eq. ((2|) is pro- 
portional to <fi 5 . Thus, we expand / around ip r — up to 
forth order in ip r , but also making / vanish at ip T = — 1 
(the condition for existence of absorbing states at ±1): 

f{A) = \{l + W (c + ofo + #r ~ Wr) > ( 4 ) 

where the real coefficients a, 6, c and d are, for conve- 
nience, defined as (primes denoting derivatives) 



c = 2/(0), o = 2/'(0)-c, 
d = /"(0)-o, 6 = -^^ + d 



(5) 



To obtain a closed equation for 4> r {t), we replace expres- 
sion (0| for f(ipk) into Eq. @ and make the substitution 
ip T = (f> r + A0 r , where we define the Laplacian operator 



Acj) r = i Y1t'/t (^ r ' — — — </>r- We then expand to 
first order in A</> r (assuming that 4> r is a smooth function 
of r in the long time limit, so that A(f> r <C <p T < 1), and 
obtain the following Langevin equation for cj) r 
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= (l-0 2 )(a0-60 3 ) 

+ [a + c + (d - 2a - 3o)0 2 ] A</> + 77, 



(6) 



with correlations for the noise 



(*')>= (?) 
|(1 - 2 )(c + d(j) 2 ) + (a - c + 2d)(t)A(t^5 r ^5(t - t'), 

where <fi denotes <fi r (t)- Equations ([6]) and ([7]) agree with 
the Langevin equation proposed in [3j , based on symme- 
try arguments, to describe order-disorder phase transi- 
tions in general models with two symmetric absorbing 
states (the noise correlation in their equation is a simpli- 
fied version of Eq.([7])). We have derived this expression 
from the microscopic dynamics, and therefore the differ- 
ent parameters of the theory have a clear interpretation 
as a function of the transition rates. The first two terms 
of Eq. can also be obtained by identifying the field 4> T 
with the average value (S T ) of the spin at site r over all 
spin configurations, and following the moments approach 
technique used by Krapivsky et al. for the original VM 
Q. This corresponds to the fl — > 00 limit, for which fluc- 
tuations are neglected and the equation for <p becomes 
deterministic. 

We now use Eq. ([6]) to gain insight into the macro- 
scopic ordering dynamics. At the mean-field (MF) level, 
where the noise term is neglected, Eq. ((6]) takes the 
form of a time-dependent Ginzburg-Landau equation [ItJ 
= DA4> - f^, with the potential V{cj>) = -f 4> 2 + 

SL ^ L <j) A — |</> 6 and D an effective diffusion constant. In 
Fig. Q] we sketch the shape of V and its associated f(ip) for 
different values of a and b. When a > (/'(0) > /(0)), 
V has two symmetric minima, thus the system coarsens 
driven by surface tension 17|. On the contrary, when 



a < 0, (/'(0) < /(0)), the minimum is at = 0, then 
the system remains in an active disordered state with 
particles continuously flipping their spins, and a global 
magnetization that fluctuates around zero. 

To illustrate our previous results, we now analyze a 
general class of 3-state models Il4j , known to exhibit 
curvature driven by surface tension, as recently shown 
in [15l |. They are composed by two external absorbing 
states S = ±1, and an intermediate state S = 0. The 
transition of a particle from S = — 1 to S = 1 happens 
in two stages. If we denote by p_,po and p+, the den- 
sities of nearest- neighboring particles in states —1,0 and 
1, respectively, the particle first switches from S = =Fl to 
S = with probability p±+po/2, and then with the same 
probability p± + po/2 from 5 = to 5 = ±1. Hence, 
disregarding the intermediate state, this model can be 
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thought as one with two absorbing states S = ±1, and 
an effective transition probability from S = — 1 to S = 1 
equals to (p + + p /2) 2 . In the same way, we believe that 
models with many intermediate states behave as equiva- 
lent 2-state models with effective transition probabilities 
that are non-linear in the local densities, and our the- 
ory can be applied. We test this by considering a 2-state 
model proposed by Abrams and Strogatz to study the 
competition between two languages [10(. The transition 
probabilities are given by — * ±1) = P±t where q 
is a positive real number, so that the q = 1 case reduces 
to the original VM, while q > 2 corresponds to mod- 
els with one or more intermediate states. In terms of 
the local field if) = 2p + — I, the flipping probability is 

f(ip) = (^~y^J ■ Calculating the coefficients defined in 
Eq. ([5]), and replacing them into Eq. ^ we obtain 



dj)_ 
~dt 



(g-1) 
3 x 1i 



[1-4?) [60+( (? -2)(< z -3)0 3 



±.[2 + {q-l){q-A)4> 2 ] A^ + to 



(8) 



When q < 1 (/'(0) < /(0)), the stable solution is <f> = 0, 
corresponding to a disordered active system. When q > I 
(/'(0) > /(0)), the stable solutions are <j> = ±1, thus the 
system orders driven by surface tension until it reaches 
one of the absorbing states ((f) = 1,-1 for all r). In 
particular, this type of ordering is observed when q = 2, 
for which the Langevin equation has a similar form as 
the one derived by Dall'Asta et al. for 3-state models 



dt 



</> 3 ) + (l-0 2 )A0 + VW 2 V, (9) 



confirming the above mentioned equivalence with a 2- 
state model with quadratic transition probability. The 
special case q = 1 corresponds to the VM, with Langevin 
equation 



dt 



Acj) + ^/l - <j> 2 - A(j) 7]. 



(10) 



Neglecting the Laplacian in the noise term, this equation 
is the same as the one suggested by Dickman et al. fl8j ) . 
Even though the potential is zero, the system still orders 
due to the presence of the Laplacian, but without surface 
tension. 

Numerical simulations of the microscopic dynamics. 
Equation ((6| shows that, at the MF level, the macro- 
scopic behavior of a particular model defined by the prob- 
ability / is only determined by the coefficients a and 6, 
that are expressed in terms of the derivatives of / around 
0k = 0. As qualitatively predicted by the MF theory, and 
confirmed in for a fix value b < b*, there is a unique 
generalized voter (GV) order-disorder transition at a crit- 
ical value acvi that separates an active stationary state 
with absolute magnetization m = for a < oqv, from a 
frozen ordered state with m = 1 for a > acv- For b > b* , 




FIG. 1: Flipping probability / (left boxes) and its associated 
potential V (right boxes) vs local field if). Curves correspond 
to coefficient values a = —0.3 (dotted), a — (solid), a = 0.3 
(dashed), for b = -0.25 (top) and b = 1.0 (bottom). For 
both values of 6, a single-well and double-well potentials are 
obtained for a < and a > respectively. 



as a is increased, a symmetry breaking Ising transition 
is observed at a value a/, followed by a directed percola- 
tion (DP) transition at a value aop > a/- For a < a/ the 
system is disordered (m = 0), whereas for ai < a < aop 
it gets partially magnetized (0 < m < 1). Above ajjp, 
the system relaxes to the fully ordered state (to = 1). 
Since we know the connection between the macroscopic 
and the microscopic dynamics expressed in Eq. ([5]) and 
Eq. (U) respectively, we now study these transitions by a 
Monte Carlo simulation of the model. This approach is 
complementary to the one followed by Hammal et al. in 
which they integrate the Langevin equation; and it also 
allows to test the field theory. 

We performed numerical simulations on a 2- 
dimensional square lattice with first-nearest neighbors 
(lst-NNs) interactions (z — 4). We used flipping prob- 
abilities that are polynomial functions of the form of 
Eq. 0, for b = —0.25, 0.5 and 3.0, and various values of 
a (see Fig. [T]). Coefficients c and d were set in order to 
make / an increasing function of ip, and to arbitrarily fix 
the point f(ip = 1) = 1. Starting from an ordered config- 
uration of down spins (initial quenching) , we flipped the 
spins of four neighboring sites at the center of the lattice 
and let the system evolve. We found that the average 
density of up spins N and the survival probability P, 
for the three values of b, decay at the critical transition 
point acv a s N ~ t v and P ~ t~ s (not shown), where 
77 ~ and 8 ~ 0.95 agree with the exponents and 1.0 
respectively, of the GV universality class [l,[l8[. Also, we 
found that the average density of interfaces has, at oqv, 
a logarithmic decay with time (p ~ 7r/(21nf)), as in the 
VM Q. Surprisely, only a GV transition was observed 
for all three values of b, in contradiction with 
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FIG. 2: GV, DP and Ising transitions on a 2D square lattice 
of side L = 400 with up to 3rd-NNs (z = 12). (a) 1/p vs time 
on a log-linear scale, for b = —0.25 and values of a around the 
GV critical point cigv — —0.1105. (b) p vs time on a log-log 
scale for b — 0.5 and values of a around the DP critical point 
a,DP — 0.2127. (c) Binder cumulant U vs a for b = 0.5. Curves 
cross at ai ~ 0.205, where U ~ 0.56, close to the universal 
value 0.6107 of the d = 2 Ising model (horizontal dashed line). 
Insets: (a) survival probability P and (b) density of up spins 
N, on a 800 2 lattice, starting from a quenched configuration. 
Curves are averages over 10 5 realizations. Dashed straight 
lines have slopes 2/-7T and —1 in (a) and its inset respectively, 
whereas the slopes are —0.45 and 0.2295 in (b). 



However, our results are in agreement with that of 
Dornic et. al. .4J where, by studying the dynamics of 
coarsening without surface tension (also with lst-NNs 
interactions), they conjectured that all models with Z^- 
symmetry and without bulk noise exhibit GV transitions. 
This apparent disagreement between theory (three tran- 
sitions) and simulations (only GV transition) comes from 
the fact that Ising dynamics is observed when there is 
bulk noise, and this happens in our model when the in- 
teraction distance is equal or larger than 2. Thus, a spin 
surrounded by 8 parallel spins can still flip if at least one 
of the four 3-rd NNs is antiparallel. Indeed, simulations 
taking up to 2-nd NNs (z = 8) also revealed GV only, but 
increasing the interactions up to 3-rd NNs (z = 12), all 
GV, Ising and DP transitions were observed. This agrees 
with the work by Droz et al. Q, in which they studied 
an absorbing Ising model in two dimensions, and found 
that extending the interaction range up to z — 12 neigh- 
bors, the transition from disorder to order splits into a 
first Ising transition that breaks the symmetry, and then 
a DP transition to the unique absorbing state selected by 
the spontaneous symmetry breaking. 

In Fig.[2]we plot the numerical results for z — 12 neigh- 
bors. We see that for b — —0.25 (Figl^a)), the decay of 
p and P at acv — —0.1105 correspond to that of a GV 
transition, whereas for b = 0.5 (Fig. [2Kb)), the transition 
to complete order happens at a value aop — 0.2127 at 



which p ~ t~ s , N ~ f and P ~ t~ s (not shown), with 
S ~ 0.45 and rj ~ 0.2295, i.e. DP critical exponents. 
In order to find the Ising transition, we calculated the 
Binder cumulants {7 = 1 — 1114/31112 (Fig. EIc)), where 
7714 and 77i2 are the fourth and second moments of the 
magnetization, as a function of a. As we can see, at the 
critical point aj ~ 0.205, the curves of the Binder cu- 
mulants for different system sizes cross each other at the 
value U ~ 0.56, similar to the universal value 0.6107 of 
the 2D Ising model. 

Conclusions. Summing up, we have derived from the 
microscopic dynamics, the Langevin equation for the 
magnetization field of general non-equilibrium spin sys- 
tems with two symmetric absorbing states. This equa- 
tion agrees with the one introduced in previous work, 
but now the dependence of the different terms on the 
flipping probability is explicitly stated. This method- 
ology allows to predict the macroscopic behavior, such 
as critical properties and ordering dynamics, by simply 
knowing the derivatives of the transition probabilities. A 
large class of models in many different disciplines can be 
studied in this way. The generalization of this approach 
to models with an arbitrary number of symmetric ab- 
sorbing states seems to be challenging. 
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